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ABSTRACT 

We present an analysis of a devised sample of Rotation Curves (RCs), aimed at 
checking the consequences of a modified f(R) gravity on galactic scales. Originally 
motivated by the the dark energy mystery, this theory may serve as a possibility of 
explaining the observed non-Keplerian profiles of galactic RCs in terms of a break- 
down of the Einstein General Relativity. We show that in general the power-law f(R) 
version could fit well the observations with reasonable values for the mass model 
parameters, encouraging further investigation on R n gravity from both observational 
and theoretical points of view. 

Key words: gravitation - dark matter - galaxies: kinematics and dynamics 



1 INTRODUCTION 

It is well-known that the RCs of spiral galaxies show a non- 
Keplerian circular velocity profile which cannot be explained 
by considering a Newtonian gravitational potential gener- 
ated by the baryonic matter (Persic, Salucci & Stel 1996). 
Current possible explanation of this controversy includes, 
among others, the postulate of a new yet not detected state 
of matter, the dark matter (Rubin 1983), a phenomeno- 
logical modification of the Newtonian dynamics (Milgrom 
1983; Brownstein & Moffat 2006; Sanders & McGaugh 2002; 
Bekenstein 2007), and higher order gravitational theories 
(originally devoted to solve the dark energy issue, see e.g., 
(Carroll et al. 2004; Capozziello et al. 2004)). 

The recent theory proposed by Capozziello, Cardone & 
Troisi 2007 (hereafter CCT), modifies the usual Newtonian 
potential generated by baryonic matter in such a way that 
the predicted galaxy kinematics and the observed one have 
a much better agreement. They consider power- law fourth 
order* theories of gravity obtained by replacing in the grav- 
ity action the Ricci scalar R with a function f(R) oc R n , 
where n is a slope parameter. The idea is that the Newto- 
nian potential generated by a point-like source gets modified 
in to 



Gm _ 1 r/ . 
■ — {1+2^ 



■1]}, 



(1) 



where f3 is a function of the slope n, and r c is a scale length 
parameter. It turns out that in this theory (3 is a univer- 



* The term comes from the fact that the generalized Einstein 
equations contain fourth order derivatives of the metric. 



sal constant while r c depends on the particular gravitating 
system being studied. In a virialized system the circular ve- 
locity is related to the derivative of the potential through 
V 2 — r d(j)(r)/dr. It is clear that (1) may help in the expla- 
nation of the circular velocity observed in spirals. 

We remark that any proposed solution to the galaxy 
RC phenomenon must not only fit well the kinematics but, 
equally important, also have best-fit values of the mass 
model parameters that are consistent with well studied 
global properties of galaxies. 

For a sample of 15 Low Surface Brightness galaxies the 
model described in CCT was fairly able to fit the RCs. How- 
ever, in our view, the relevance of their finding is limited by 
the following considerations: 

• the sample contains several objects whose RCs are not 
smooth, symmetric and extended to large radii 

• the sample contains only Low Surface Brightness galax- 
ies while a wider sample is desirable 

• the universal parameter n is not estimated by the anal- 
ysis itself but it is taken from other observations. 

In the present work we generalize the results of 
Capozziello et al. 2007 and test a wider and fairer sample of 
spirals, improving the analysis methodology. Our goal is to 
perform a check of their model on galactic scales in order to 
investigate its consistency and universality. 

The plan of this article is the following: in Sect. 2 we 
briefly summarize the main theoretical results described 
in CCT relevant for the analysis of our sample. In Sect. 3 
we present our sample and methodology of analysis. In 
Sect. 4 the results are presented and finally the conclusions 
in Sect. 5. 
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2 NEWTONIAN LIMIT OF F(R) GRAVITY 

The theory proposed by Capozziello et al. 2007 is an example 
of f(R) theory of gravity (Nojiri & Odintsov 2007; Carloni et 
al. 2005) . In these theories the gravitational action is defined 
to be: 



-g [f(R)+Cr, 



(2) 



where g is the metric determinant, R is the Ricci scalar and 
C m is the matter Lagrangian. They consider: 



f(R) = f R n 



(3) 



where fo is a constant to give correct dimensions to the 
action and n is the slope parameter. The modified Einstein 
equation is obtained by varying the action with respect to 
the metric components. 

Solving the vacuum field equations for a Schwarzschild- 
like metric in the Newtonian limit of weak gravitational 
fields and low velocities, the modified gravitational poten- 
tial for the case of a point-like source of mass m, is given by 
(1), where the relation between the slope parameter n and 
(3 (see detailed calculation in CCT) is given by: 



Yin 2 



7n-l- V36n 4 + 12n 3 - 83n 2 + 50n + 1 



(4) 



6n 2 - 4n + 2 

Note that for n — 1 the usual Newtonian potential is recov- 
ered. The large and small scale behavior of the total poten- 
tial constrain the parameter (3 to be < (3 < 1. 

The solution (1) can be generalized to extended systems 
with a given density distribution p(r) by simply writing: 



4>{r) = 



-G 



|r — r 

4>n(v) +0c(r), 



-1]} 



(5) 



where 4>n(v) represents the usual Newtonian potential and 
4>c(r) the additional correction. In this way, the Newtonian 
potential can be recuperated when (3 = 0. The solution for 
the specific density distribution relevant for spiral galaxies 
is described in the following paragraph. 



3 DATA AND METHODOLOGY OF THE TEST 

We selected two samples of galaxies: a first with 15 galaxies, 
called Sample A, that represents the best available RCs to 
study the mass distribution of luminous and/or dark matter, 
and it has been used in works concerning modifications of 
gravity and the core/cusp controversy (Corbelli & Salucci 
2007; Gentile et al. 2004; Frigerio Martins & Salucci 2007). 

This sample includes nearby galaxies of different Sur- 
face Brightness: DDO 47 (Gentile et al. 2005); ESO 116-G12, 
ESO 287-G13, NGC 7339, NGC 1090 (Gentile et al. 2004); 
UGC 8017, UGC 10981, UGC 11455 (Vogt et al. 2004); M 
31, M 33 (Corbelli & Salucci 2007); IC 2574 (Martimbeau, 
Carignan & Roy 1994), NGC 5585 (Cote, Carignan & San- 
cisi 1991), NGC 6503 (Wevers, van der Kruit & Allen 1986), 
NGC 2403 (Fraternali et al. 2002), NGC 55 (Puche, Carig- 
nan & Wainscoat 1991). This sample is the most suitable for 
a fair test of theories like the one of Capozziello et al. 2007: 

• The RCs are smooth, symmetric and extended to large 
radii. 



• The galaxies present a very small bulge so that it can 
be neglected in the mass model to a good approximation. 

• The luminosity profile is well measured and presents a 
smooth behavior 

• The data are uniform in quality up to the maximal radii 
of each galaxy. 

Let us notice that in some of these galaxies H a and HI RCs 
are both available and in these cases they agree well where 
they coexist. 

We also considered a second sample called Sample B 
consisting of 15 selected objects from Sanders & McGaugh 
2002 that has been used to test MOND. This sample con- 
sists of the following galaxies: UGC 6399, UGC 6983, UGC 
6917, NGC 3972, NGC 4085, NGC 4183, NGC 3917, NGC 
3949, NGC 4217, NGC 3877, NGC 4157, NGC 3953, NGC 
4100 (Tully et al. 1996; Verheijen & Sancisi 2001); NGC 
300 (Puche, Carignan & Bosma 1990); UGC 128 (van der 
Hulst et al. 1993). Although these galaxies do not fulfill 
all the requirements of Sample A we have analyzed them 
for completeness sake. The properties of the galaxies of the 
two samples are listed in Table 1. Notice that the theory of 
Capozziello et al. 2007 requires an analysis with a sample of 
high quality galaxies, as described above, where each lumi- 
nous profile plays an important role, whereas this is not the 
case in MOND. 

We decompose the total circular velocity into stellar 
and gaseous contributions. Available photometry and radio 
observations show that the stars and the gas in our sample 
of galaxies are distributed in an infinitesimal thin and circu- 
lar symmetric disk. While the HI surface luminosity density 
distribution Yi gas {r) gives a direct measurement of the gas 
mass, optical observations show that the stars have an ex- 
ponential distribution: 



Eu(r) = (M D /27iR 2 D )e- r/RD , 



(6) 



where Md is the disk mass and Rd is the scale length, the 
latter being measured directly from the optical observations, 
while Md is kept as a free parameter of our analysis. 

The distribution of the luminous matter in spiral galax- 
ies has to a good extend cylindrical symmetry, hence using 
cylindrical coordinates, the potential (5) reads 



(r) 



E(r ; ) is the surface density distribution of the stars, given by 
(6) , or of the gas, given by an interpolation of the HI data 
points up to the last measured point. f3 and r c are free pa- 
rameters of the theory, with the latter galaxy dependent. We 
neglected the gas contribution to the mass density for radii 
larger than the last measured point, however we checked the 
goodness of this approximation by extending the distribu- 
tion with a different kind of decreasing smooth curves and 
realized that error made in the truncated approximation is 
small enough to be neglected. 

Denning k 2 = , 4r r , - , we can express the distance be- 

(r+r ) z 

tween two points in cylindrical coordinates as |r — r'| = 
(r + r) 2 (l — k 2 cos 2 (0/2)). The derivation of the circular ve- 
locity due to the marked term of equation (7), that we call 
(j)f3 (r), is now direct: 
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r -r Mr) = -2 /3 " 3 r c - /3 n a {p - 1) G I(r), (8) 
ar 

where the integral is denned as 



(9) 



X(r)= [ drW'^-—^k 3 - (3 T i {r')T{r), 
Jo 2 

with T(r) written in terms of confluent hyper-geometric 
function: T(r) = 2(r + r') 2 Fi[|, 1, /c 2 ] + [(k 2 - 2)r' + 

The total circular velocity is the sum of each squared 
contribution: 

VCCT(T) — Vn, stars + VN,gas + V<3, stars + Vc,gas (10) 

where the stars and gas subscripts refer to the different 
contributions of luminous matter to the total potential (5). 
The N and C subscripts refer to the Newtonian and the 
additional correction potentials. 
Let us recall that we can write 



s (r) = (GM D /2R D )x 2 B(x/2), 



(11) 



where x = v/Rd, G is the gravitational constant and the 
quantity B = IqKq — I\K\ is a combination of Bessel func- 
tions (Freeman 1970). 

Galaxies UGC 8017, M 31, UGC 11455 and UGC 10981 
presents a very small amount of gas and for this reason it has 
been neglected in the analysis. Notice that the correction to 
the Newtonian potential in equation (1) may be negative 
and this would lead to a negative value of Vq. In Figures 1 
and 2 however the velocities Vc are shown only in the ranges 
of r where their square are positive. 

In a first step, the RCs are x 2 best-fitted with the fol- 
lowing free parameters: the slope (f3) and the scale length 
(r c ) of the theory, and the gas mass fraction (f ga s) related 
to the disk mass simply by M D = M gas (l - fgas)/ fgas- The 
errors for the best fit values of the free parameters are calcu- 
lated at one standard deviation with the xled + 1 m l e - From 
the results of these fits we get a mean value of (3 = 0.7 ±0.25 
(n ~ 2.2). In the second step we redo the best-fit fixing the 
slope parameter at f3 — 0.7 keeping as free parameters only 
r c and f gas - Notice that in a previous paper (Capozziello et 
al. 2006), a mean value of (3 = 0.58 ± 0.15 (n ~ 1.7) has 
been obtained, perfectly compatible with our result. This 
parameter however, is well constrained from SNela observa- 
tions to be (3 — 0.87 (n c± 3.5), also compatible with our 
measurements. In our analysis the value j3 = 0.7 is the most 
favorable for explaining the RCs: different values of (3 from 
the one we adopt here lead to worse performance. 



• not clear comprehension for the big variation of values 
for the scale length parameter (0.005 kpc< r c <1.53 kpc). 

The residuals of the measurements with respect to the 
best-fit mass model are in most of the cases compatible with 
the error-bars, see Figures 1 and 2, though three galaxies 
show significant deviations: NGC 6503, NGC 2403 and M 
33. We also find acceptable values for the B-band mass- 
to-light ratio parameter for most of the galaxies, for which 
we should have approximately 0.5 < Tf < 6 and a positive 
correlation between B- luminosity t and Tf (Salucci & Persic 
1999): 

M d (Lb) ~ 3.7 xlO 10 

x [(^) 1 - 23 ^(L B )+O.O95(^)°- 98 ]M , (12) 

where L 10 = 10 10 L B q and g(L B ) = exp[-0.87 x (logj^ - 
0.64) 2 ]. In detail we find discrepancies for NGC 55, UGC 
8017, NGC 3972, NGC 4085 and NGC 4183. Values for the 
scale length parameter (r c ) are in general smaller for less 
massive galaxies and bigger for more massive ones. We ob- 
tained a Newtonian fit for UGC 10981, as shown by the 
exceedingly large value for r c , see Figure 1. 

The model analyzed in this article yields better results 
on galactic scales than A Cold Dark Matter models, where in 
the latter these galaxies have serious problems like marginal 
fits and unreasonable values for the stellar mass-to-light ra- 
tio, see e.g., Frigerio Martins & Salucci 2007 and Gentile et 
al 2004. 



5 CONCLUSIONS 

We have investigated the possibility of fitting the RCs of 
spirals with a power-low fourth order theory of gravity of 
Capozziello et al 2007, without the need of dark matter. 
We remark the relevance of our sample that contains ob- 
jects in a large range of luminosity and with very accurate 
and proper kinematic. We find in general a reasonable agree- 
ment, with some discrepancies, between the RCs and the 
Capozziello et al. 2007 circular velocity model, encouraging 
further investigations from the theoretical point of view. 

Acknowledgments. We warmly thank S. Capozziello, J. 
Miller and T. Sotiriou for useful discussions. This research 
was supported by CAPES-Brasil (C.F.M). 



4 RESULTS 

We summarize the results of our analysis in Figures 1 and 
2 and Table l"L In general we find for all galaxies: 

• the velocity model Vcct well fitting the RCs 

• acceptable values for the stellar mass-to-light ratio 

• too vast range for values of the gas fraction (0% < f g < 
100%) 



t Numerical codes and data used to obtain these results can be 
found at the address http://people.sissa.it/~martins/home.html 



t = Md/Lb] Md is the disk mass and Lb is the B-band 
galaxy luminosity 
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Figure 1. Sample A: The solid line represents the best-fit total circular velocity Vqct- The dashed and dotted lines are the Newto- 
nian contributions from the gas and the stars, while the dot-dashed represents their sum. The long-dashed line is the non-Newtonian 
contribution of the gas and the stars to the model. Below the RCs, we plot the residuals iy o bs — Vcct)- See Table 1 for details. 
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Figure 2. Sample B: Best-fit curves superimposed to the data from selected objects from Sanders & McGaugh 2002. See Figure 1 for 
details. 
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Table 1. Properties and parameters of the mass model of the analyzed Samples ((3 = 0.7). From left to right, the columns read: name 
of the galaxy, Hubble type as reported in the NED database, adopted distance in Mpc, B-band luminosity in 10 9 L#q, disk scale length 
in kpc, gas mass in 1O 9 M0 until last measured point, gas fraction in %, disk mass in 1O 9 M0, scale length CCT parameter in kpc, 
mass-to-light ratio in T5, and X^ ed - The galaxies are ordered from top to bottom with increasing luminosity. 



Galaxy 


Type 


D 


L B 


Rd 


M gas 


fgas 


M D 


r c 


1 ★ 


Y 2 












Sample 


A 










DDO 47 


IB 


4 


0.1 


0.5 


2.2 


96±1 


0.01 


0.005 


0.1 


0.5 


IC 2574 


SABm 


3 


0.8 


1.78 


0.52 


79±12 


0.14 


0.017±0.003 


0.2 


0.8 


NGC 5585 


SABc 


6.2 


1.5 


1.26 


1.45 


58±3 


1 


0.038±0.004 


0.7 


1.4 


NGC 55 


SBm 


1.6 


4 


1.6 


1.3 


84±7 


0.24 


0.024±0.004 


0.06 


0.14 


ESO 116-G12 


SBcd 


15.3 


4.6 


1.7 


21 


50 


2.1 


0.05±0.01 


0.5 


1.2 


NGC 6503 


Sc 


6 


5 


1.74 


2.3 


18±0.7 


10.6 


0.21±0.014 


2.1 


18 


M 33 


Sc 


0.84 


5.7 


1.4 


3.7 


53±2 


3.3 


0.075±0.004 


0.58 


25 


NGC 7339 


SABb 


17.8 


7.3 


1.5 


6.2 


2.8±0.2 


22 


0.41±0.07 


3 


2.3 


NGC 2403 


Sc 


3.25 


8 


2.08 


4.46 


27±0.9 


12.1 


0.21±0.015 


1.5 


19 


M 31 


Sb 


0.78 


20 


4.5 


- 


- 


180±70 


1.53±0.19 


9 


3.4 


ESO 287-G13 


Sbc 


35.6 


30 


3.3 


14 


25±1 


41 


0.48±0.05 


1.4 


3.2 


NGC 1090 


Sbc 


36.4 


38 


3.4 


100 


18±1 


47 


0.59±0.04 


1.2 


0.9 


UGC 8017 


Sab 


102.7 


40 


2.1 


- 


- 


9.1±0.3 


0.01±0.01 


0.23 


5.2 


UGC 11455 


Sc 


75.4 


45 


5.3 


- 


- 


74±3 


0.14±0.01 


1.6 


5 


UGC 10981 


Sbc 


155 


120 


5.4 


- 


- 


460±200 


-10 11 


3.8 


4.9 












Sample 


B 










UGC 6399 


Sm 


18.6 


1.6 


2.4 


1 


23±3 


3.3 


0.1±0.03 


2 


0.1 


NGC 300 


Scd 


1.9 


2.3 


1.7 


1.3 


39±4 


2 


0.052±0.010 


0.87 


0.43 


UGC 6983 


SBcd 


18.6 


4.2 


2.7 


4.1 


24±2 


13 


0.46±0.1 


3.1 


0.88 


UGC 6917 


SBd 


18.6 


4.4 


2.9 


2.6 


U±l 


16 


0.71±0.17 


3.6 


0.47 


UGC 128 


Sd 


60 


5.2 


6.4 


10.7 


32±5 


23 


0.39±0.11 


4.4 


0.1 


NGC 3972 


Sbc 


18.6 


6.7 


2 


1.5 


39±3 


2.5 


0.025±0.004 


0.37 


0.1 


NGC 4085 


Sc 


18.6 


6.9 


1.6 


1.3 


44±4 


1.7 


0.014±0.003 


0.25 


1 


NGC 4183 


Scd 


18.6 


9.5 


1.4 


4.9 


60±6 


3.2 


0.09±0.023 


0.3 


0.33 


NGC 3917 


Scd 


18.6 


11 


3.1 


2.6 


22±1.5 


9.2±0.9 


0.098±0.014 


0.8 


1 


NGC 3949 


Sbc 


18.6 


19 


1.7 


4.1 


19±2.2 


17 


0.22±0.06 


0.9 


0.25 


NGC 4217 


Sb 


18.6 


21 


2.9 


3.3 


6.1±0.7 


52 


0.55±0.15 


2.5 


0.38 


NGC 4100 


Sbc 


18.6 


25 


2.5 


4.4 


13±1.5 


28 


0.20±0.03 


1.1 


1.52 


NGC 3877 


Sc 


18.6 


27 


2.8 


1.9 


7.3±0.8 


24 


0.2±0.04 


0.9 


0.75 


NGC 4157 


Sb 


18.6 


30 


2.6 


12 


26±2.6 


33 


0.25±0.04 


1.1 


0.53 


NGC 3953 


SBbc 


18.6 


41 


3.8 


4 


2.8±0.18 


140 


1.9±0.5 


3.4 


0.78 
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